Photonic realization of nonlocal memory effects and non-Markovian quantum probes 
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The study of open quantum systems is important for fundamental issues of quantum physics 
as well as for technological applications such as quantum information processing. Recent develop- 
ments in this field have increased our basic understanding on how non-Markovian effects influence 
the dynamics of an open quantum system, paving the way to exploit memory effects for various 
quantum control tasks. Most often, the environment of an open system is thought to act as a sink 
for the system information. However, here we demonstrate experimentally that a photonic open 
system can exploit the information initially held by its environment. Correlations in the environ- 
mental degrees of freedom induce nonlocal memory effects where the bipartite open system displays, 
counterintuitively, local Markovian and global non-Markovian character. Our results also provide 
novel methods to protect and distribute entanglement, and to experimentally quantify correlations 
in photonic environments. 

PACS numbers: 03.65.Yz, 42.50.-p, 03.67.-a 



Realistic quantum systems interact and exchange in- 
formation with their surroundings. The engineering of 
the decoherence and the flow of information between an 
open quantum system and its environment has recently 
allowed, e.g., to drive quantum computation by dissipa- 
tion [1 , to control entanglement and quantum phases 
in many-body systems [2HI], to create an open system 
quantum simulator [5 , and to control the transition from 
Markovian dynamics to a regime with quantum mem- 
ory effects [6 j. Moreover, there has been recently rapid 
progress in understanding various fundamental issues on 
quantum non-Markovianity [7UT3]. 

In general, the power of quantum information process- 
ing is based on the presence of quantum correlations in 
multipartite systems. However, as any realistic quantum 
system, those systems are vulnerable to the detrimental 
effects of the environment such as the loss of coherence 
and information. Depending on its physical composition, 
an open multipartite system can interact either with a 
single common environment or each of the sub-parties 
can interact locally with their own reservoirs. The lat- 
ter case, which is in the focus of our paper, is a common 
scenario, e.g., in quantum information processing and en- 
tanglement dynamics [14], and in energy transport in bi- 
ological organisms [T5] . 

The central question to be studied here is, can a com- 
posite open quantum system exploit during its time evo- 
lution the information carried initially by the environ- 
mental degrees of freedom, i.e., the information to which 
it does not have access initially? On the basis of a recent 
theoretical prediction [16 , we demonstrate experimen- 
tally that initial correlations between local parts of the 
environment lead to nonlocal memory effects and allow 
to protect and maintain genuine quantum properties of 



an open system. We will show further that the non- 
Markovian dynamics of photonic open systems provides 
a diagnostic tool for characteristic features of their envi- 
ronments. 

The physical system and experimental setup. Our ex- 
perimental open system consists of the polarization de- 
grees of freedom of a pair of entangled photons. The 
photon pair is created by spontaneous parametric down- 
conversion after which the photons travel along different 
arms i = 1,2 and move through different quartz plates 
of variable thickness (see Fig. [T]). The quartz plates act 
as birefringent media which lead to a local coupling be- 
tween the polarization degrees of freedom of the photon 
in arm i and its frequency (mode) degrees of freedom 
forming its local environment. This local interaction can 
be described by a local unitary operator which is defined 
by Ui(U)\X) <8> \<Ji) = e in ^ u \\) <g> |^), where |A) <g> \oOi) 
denotes the state of the photon in arm i with polariza- 
tion A = H, V (horizontal or vertical) and frequency ooi. 
The refraction index for photons with polarization A is 
denoted by n\, and U represents the interaction time. 

Any pure initial state for the polarization degrees of 
freedom of the photon pair, which forms our bipartite 
open system, can be written as 1^12) = a \ HH)+b \ HV) + 
c \ VH) + d \VV). The total system initial states are as- 
sumed to be product states 

|*(0)) = |^i2> ® J doo x doo 2 g(oo x ,oo 2 ) ^1,^2) , 

where g(ooi : oo 2 ) 1S the probability amplitude for the pho- 
ton in arm 1 to have frequency oo\ and for the photon in 
arm 2 to have frequency oo 2 ? with the corresponding joint 
probability distribution P (001,002) — \g(w\, ^2) | 2 • 
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or by tracing over photon 1 to find the state of photon 2: 
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FIG. 1: (color online). Experimental setup. Maximally en- 
tangled photon pairs are generated by parametric down con- 
version (PDC) with an ultraviolet (UV) pump pulse whose 
spectral width can be controlled. Then a series of quartz 
plates are added in arm 1 and arm 2 to realize the local de- 
phasing channels. The final two-photon polarization state is 
analyzed by state tomography. Key to the components: NC - 
nonlinear crystal, FP - Fabry-Perot cavity, DM - dispersion 
medium, HWP - half wave plate, QWP - quarter wave plate, 
IF - interference filter, QP - quartz plate, PBS - polarizing 
beamsplitter, and SPD - single photon detector. 



Open system dynamics. The time evolution of the 
open system can be represented by means of the dy- 
namical map &\ 2 El which maps the initial polariza- 
tion state pi2(0) = | '012) ('0i2| to the polarization state 
state pi 2 (t) = $12(^12(0)) at time t. Calculating the to- 
tal system time evolution and taking the trace over the 
environment yields [16] 
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The dynamics represents a pure dephasing process whose 
decoherence functions can be expressed completely in 
terms of the Fourier transform of the joint probability 
distribution 



G(h,t 2 ) = J duJiduj 2 P(uJi, 



cu 2 )e 



-iAn(uj\t-L+0J2t2) 



(2) 



as = G(*i,0), K 2 (t) = G(0,t 2 ), «i 2 (*) = G(t u t 2 ) 

and Ai 2 (t) = G(ti, — £2), where An = riy — nn denotes 
the birefringence. The function K>i(t) describes the de- 
coherence of the polarization state of photon i which is 
induced by the local coupling to its frequency environ- 
ment caused by the quartz plate in arm i. This can be 
seen by taking the trace over the polarization state of 
photon 2 to obtain the polarization state of photon 1: 
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Nonlocal dynamical map. The nonlocal character of 
the dynamical map Q\ 2 is controlled by the decoherence 
functions K\ 2 (t) and A 12 (t) in Eq. 0. As a matter of 
fact, the map &\ 2 is a product of local dynamical maps, 
i.e., $[ 2 = §\ (g) $2 ^ an d only if K 12 (t) = Ki(t)n 2 (t) 
and Ai2(t) = ni{t)K^{t). In particular, this implies that 
frequency correlations in the initial state do not influ- 
ence the local dynamics of each photon and cannot be 
detected by observing the local polarization states. How- 
ever, the evolution of the composite two-photon polariza- 
tion state is determined by the frequency correlations in 
the environment which lead to memory effects and global 
non-Mar kovian dynamics. This also means that, coun- 
terintuitively, by adding degrees of freedom to an open 
system one can change its dynamics from the Markovian 
to the non-Mar kovian regime. 

Trace distance dynamics. To quantify non- 
Markovianity we use the trace distance based 
measure [11] which is defined by M($) = 
max PA,B(o) f a> o dt(7 ( t ^PA,B(0)). Here, a(t, Pa,b(0)) = 
■^D(pA(t)j pB(t)) and the trace distance D between two 
states p A and p B is given by D(p A ,p B ) = ^ti\p A - Pb\- 
The amount of non-Markovianity is therefore equal to 
the total increase of the trace distance during the time 
evolution and quantifies the total information flow from 
the environment to the system. This measure has been 
used recently in several theoretical and experimental 
contexts, see, e.g. Refs. [6| fT5| fT8| fT9] . 

For anticorrelated frequencies uji the pair of the Bell- 
states \i/>f 2 ) = i\HH) ± \ VV)) /y/2 maximizes the in- 
crease of the trace distance and thus determines the non- 
Markovianity measure M [16 ]. These states are created 
by using spontaneous parametric down conversion as is 
illustrated in Fig. [l] A femtosecond pulse (the duration 
is about 150 fs and the operation wavelength is at 780 nm, 
with a repetition rate of about 76 MHz) generated from 
a Tiisapphire laser is frequency doubled to pump two 
lmm-thick beamlike cut beta barium borate (BBO) crys- 
tals [20] creating the two-photon entangled state. With 
the help of QWPl we can easily tune the entangled state 
between |^+) (QWPl is set at 0°) and \^ 2 ) (QWPl 
is set at 90°). The size of the anticorrelations between 
the photon frequencies is controlled in the experiment 
by the spectrum of the UV pump pulse. Since energy is 
conserved in the down conversion process, decreasing the 
spectral width of the pump pulse decreases the uncer- 
tainty of the sum of the frequencies of the photons, and 
hence increases the anticorrelation between the frequen- 
cies. We use four different pulse widths, see Supplemental 
Material at [URL will be inserted by publisher] for more 
details. It is also worth noting that even though our 
experimental setup is based on downconversion, which 
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FIG. 3: (color online). Trace distance dynamics of the open 
system with different dispersions of the pump pulse. The 
quartz plate on the path of the pump pulse has thickness (a) 
(b) 13.19 mm (c) 39.57 mm. In all cases K — —0.55 and the 
measured values of non-Mar kovianity are J\f — 0.14, 0.13, 0.13. 
The trace distance dynamics remains the same irrespective of 
the character of the correlations. See Supplemental Material 
at [URL will be inserted by publisher] for further details. 



FIG. 2: (color online). Trace distance dynamics for different 
values of the spectral width of the pump pulse. The widths 
and the amounts of non-Markovianity are: (a) 5 = 0.18 nm, 
Af = 0.48, (b) 6 = 0.52 nm, N = 0.23, (c) 6 = 0.73 nm, 
AT = 0.14, (d) 6 = 1.89 nm, J\f = 0.02. The x-axis represents 
time measured by the total effective path difference between 
the horizontal and vertical photons caused by quartz plates in 
arms 1 and 2 (Ao = 780 nm). The y-axis is the trace distance 
between the time-evolved initial open system states \ipi 2 ) and 
\ipi 2 )- We first add quartz plates in arm 1. When the to- 
tal effective path difference between horizontal and vertical 
photons equals 199Ao, we add quartz plates in arm 2. The 
solid line shows the fit using the theoretical result of Eq. (J3| . 
The error bars are due to the counting statistics. See Sup- 
plemental Material at [URL will be inserted by publisher] for 
more details on the counting statistics and the fitting of the 
theoretical curves. 



is a commonly used tool to prepare specific two-photon 
states in quantum optical experiments, the physical phe- 
nomena that we describe and realize experimentally was 
theoretically predicted only very recently [16]. 

Assuming a Gaussian joint frequency distribution 
P{(jJi,uj2) with identical single frequency variance C and 
correlation coefficient K, the time evolution of the trace 
distance corresponding to the optimal Bell-state pair is 
found to be [16] 



D(t) exp 



1 



An 2 C(^+^-2|if|M 2 ) 



(3) 



For uncorrelated photon frequencies, i.e. for uncorrelated 
local environments we have K = and the trace dis- 
tance decreases monotonically, corresponding to Marko- 
vian dynamics. However, as soon as the frequencies 
are anticorrelated, K < 0, the trace distance is non- 
monotonic which signifies quantum memory effects and 
non-Markovian behavior. 

Consecutive interactions with local environments. In 
our first experiment the local decoherence induced by 
the quartz plates in arm 1 and arm 2 act consecutively, 
and we change the size of the initial correlations between 



the local reservoirs by tuning the spectral width of the 
pump pulse. The results for the open system's trace dis- 
tance dynamics are displayed in Fig. [2] which shows how 
the initial environmental correlations induce and influ- 
ence the quantum non-Markovianity. First, when only 
the quartz plate in arm 1 is active, the trace distance 
decreases monotonically demonstrating that information 
flows continuously from the system to the environment. 
When subsequently the quartz plate in arm 1 becomes 
inactive and the quartz plate in arm 2 active, the trace 
distance increases again highlighting non-Markovian be- 
havior and a reversed flow of information from the en- 
vironment back to the open system. It is important to 
note that this increase of the trace distance implies a 
revival of the entanglement between the photon polar- 
izations. Whilst there exists earlier literature, e.g., on 
how the entanglement is transferred from the open sys- 
tem to the environment [21] , or on entanglement dynam- 
ics when the qubits interact with local non-Markovian 
reservoirs [14], our results are fundamentally different. 
We demonstrate that the entanglement can revive due 
to the nonlocal character of the dynamical map despite 
of the fact that the local interaction of the qubits with 
their environments is completely Markovian. A further 
important aspect of the experimental scheme is that it 
enables to measure the frequency correlation coefficient 
K of the photon pairs by fitting the theoretical predic- 
tion of Eq. ([3| to the experimental data. The fits shown 
in Fig. [2] corresponding to panels (a)-(d) yield the val- 
ues K = -0.92, -0.66, -0.55, and -0.17. Thus we see 
that the open system (polarization degrees of freedom) 
functions as a quantum probe which allows us to gain 
nontrivial information on the correlations in the envi- 
ronment (frequency degrees of freedom). To the best of 
our knowledge, our results present the first quantitative 
measurement of the correlation coefficient without recon- 
structing the joint frequency distribution. 

Classical vs. quantum initial correlations between the 
local environments. The open system dynamics of Eq. ([!]) 
depends solely on the decoherence functions ^i, ^2, ^12 
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FIG. 4: (color online). Trace distance dynamics for five 
different quartz plate configurations and a fixed correlation 
coefficient of K — —0.92. For all curves the maximal length 
of the quartz plates in both arms is 199Ao, giving a maximal 
total effective path difference of 398Ao- From top to bottom, 
the curves show the results for the cases, when we start to add 
quartz plates to arm 2 when the length of the plate in arm 
1 has reached the length of 0, 75, 100, 150, 199Ao respectively, 
see Supplemental Material at [URL will be inserted by pub- 
lisher] for details. Hence, in the top curve the quartz plates 
are added simultaneously to both arms while in the bottom 
curve they act consecutively. The results demonstrate that 
the final state is the same while in terms of the information 
flow the behaviour is completely different. For the top (red) 
curve information flows slowly out of the system with a small 
rate and, in fact, the decoherence would completely halt in the 
ideal ideal case of fully anticorrelated frequencies (K = —1.0). 
In the other extreme, for the bottom (violet) curve informa- 
tion first flows out of the system with a high rate and later 
on returns back nearly completely indicating strong memory 
effects. The solid lines are fits on the basis of Eq. see 
Supplemental Material at [URL will be inserted by publisher] 
for more details. 



and A12 which in turn are completely determined by the 
Fourier transform ([2| of the joint frequency distribution 
P(lji,lj2)- Therefore, also the dynamical map of the 
open system depends only on P{lji,lj2) and not on the 
probability amplitudes g(co>i, 0J2). This fact leads to the 
conclusion that the nonlocal memory effects, and the re- 
vival of entanglement observed above, do not depend on 
the specific phase relations between different frequency 
components, and consequently, whether the initial fre- 
quency correlations between the photons are of quantum 
or classical character. To test experimentally whether 
the nature of the correlations plays a role in inducing 
the memory effects, we insert quartz plates also on the 
path of the pump pulse, thereby changing the dispersion 
relations of the pump pulse while keeping fixed the cor- 
relation coefficient K. The results in Fig. [3] show that 
the trace distance dynamics remains the same in this 
case and only depends on the correlation coefficient K of 



the photon frequencies. Thus we see that nonlocal mem- 
ory effects and the revival of the entanglement between 
the polarization states of the photons can be induced by 
purely classical correlations between local environments 
of the open system. 

Different interaction configurations. In the previous 
experiments the local dephasing environments acted con- 
secutively. In our final experimental series we reposition 
the quartz plates in arm 2 so that eventually the dephas- 
ing is active simultaneously in both arms (see Supple- 
mental Material at [URL will be inserted by publisher] 
for the schematic view of the configurations). The results 
are shown in Fig. [4j We observe that the final state of 
the open system is always the same when the total ef- 
fective path difference is equal to 398 Ao- However, the 
amount of non-Markovianity and the flow of information 
from the environment to the open system decreases with 
increasing overlap of the action of the quartz plates. This 
demonstrates that initial correlations between the local 
environments and the reordering of the quartz plates al- 
low to engineer the information flow between local sys- 
tems while keeping fixed the final state. It is also impor- 
tant to note that when the local dephasing channels act 
simultaneously on our open bipartite system (the top red 
curve in Fig.[4|, the initial correlations between the com- 
posite environments protect the entanglement between 
the qubits. One can see this also from Eq. (|3| which 
shows that for t\ = £2 and K = — 1 the trace distance 
remains constant and that there is thus no outflow of in- 
formation from the system for the considered Bell-states. 

Conclusions. We have realized nonlocal quantum dy- 
namical maps of a photonic system and demonstrated 
experimentally how initial quantum or even classical cor- 
relations between the local environments of a composite 
open system induce nonlocal memory effects. While the 
global dynamics of the open system is non-Markovian, its 
local subsystems behave perfectly Markovian. The obser- 
vation of the open system dynamics yields also important 
information about the environment. In fact, we have 
seen that a measurement of the non-Markovian evolu- 
tion of the two-photon polarization state leads to a novel 
method for the determination of the correlation coeffi- 
cient of the photon frequencies. Thus, the open system 
dynamics can be viewed as a non-Markovian quantum 
probe which enables to extract nontrivial information 
on characteristic features of the environment. This fact 
as well as our demonstration that initial correlations in 
the environment can protect the entanglement within an 
open system, could find important applications in quan- 
tum technologies and control strategies, and in the devel- 
opment of novel diagnostics tools for complex quantum 
systems. 
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Controlling the widths of the pump pulse 
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The width of the pump spectrum controls the fre- 
quency anticorrelations between the pair of photons cre- 
ated in the down conversion process. By using filters and 
fused silica plates, we can control the width of the pump 
pulses which are depicted in Fig. [5) The laser source is 
filtered to 3nm (FWHF) and then passes through the 
frequency doubler (1.5mm thick BiBaC^). Hence, the 
bandwidth of the UV pulse is about 0.52 nm as shown 
in Fig. [5] (b). To obtain a sharper spectrum, as shown 
in Fig. pRa), we insert a thin fused silica plate, which is 
0.05 mm thick and coated with a partial reflecting coat- 
ing on each side, with approximately 75% at 390 nm. For 
the spectra displayed in Fig. [5] (c) and (d), we have no 
filter before the doubler and no fused silica plate after the 
doubler, and for Fig.[5](d) we also replace the doubler by 
a 0.3mm-thick BBO. 



Counting statistics 

The details of the counting statistics and pump power 
for the experimental results of Fig. 2 of the main text 
are as follows. For Fig. 2 (a) the pump power is about 
2.4 mW and the total coincident count rate on the basis 
vectors HH, HV, VH, VV is about 18, 000 in 10 seconds, 
(b) 10 mW and the rate is about 35, 000 in 4 seconds, (c) 
20 mW and the rate is about 35,000 in 2 seconds, and 
(d) 20 mW and rate is about 36, 000 in 4 seconds. Please 
note that the coincident count rate always depends also 
on the bandwidth of UV pump source. 



Fitting the experimental data 



FIG. 5: (color online). The frequency spectra of the ultra- 
violet (UV) pump pulses used for the down conversion pro- 
cess. Full widths at half maximum (FWHM) of the pulse 
spectra are (a) 5 = 0.18 nm, (c) S = 0.52 nm, (c) S = 0.73 nm, 
and (d) 5 — 1.89 nm. The solid lines represent Gaussian fits 
of the experimental data. 



process, we can find A and B which agree with the exper- 
imental data. Then in the second part, when the quartz 
plate in arm 2 is active, we use the function g(x) = 
Aexp [-B (199 2 + (x - 199 2 ) - 2\K\199(x - 199))] . 
Here, 199 corresponds to the effective path difference 
after the first quartz plate. Since A and B are deter- 
mined by the first fit, we can easily determine K from 
the second fit. 

For Fig. 4, again keeping in mind Eq. (3), we can do 
the fitting in the following way. First, we use the func- 
tions f(x) and g{x) mentioned above to make a fit of 
the lowest curve and determine Ai, B\, and K\. Then 
we use a function G(x) = Aexp [-B (x 2 /2 - \K\x 2 /2)] 
to fit the top curve of the figure. From here, we de- 
termine A2, B<i, and K2. Then we determine Ao = 
(A 1 + A 2 )/2, B = (B 1 + B 2 )/2, and K = [K x + K 2 )/2 
which are then used to plot the functions F(ti,t 2 ) = 
A exp [-B {t\ + tl - 2\Ko\ht 2 )] in the figure. 



Trace distance dynamics with different dispersions 
of the pump pulse 



Employing Eq. (3), the theoretical fits of the experi- 
mental data presented in Figs. 2 and 3 of the main text 
have been carried out as follows. We first use the func- 
tion f{x) = Aexp (—Bx 2 ) to make a fit for the process 
when the quartz plates are added to arm 1. From this 



For the experimental results of Fig. 3 of the main text, 
the FWHM of the pump pulse is about 0.73 nm, the 
pump power is about 20 mW, and the total coincident 
count rate on the basis vectors HH, HV, VH, VV is about 
35, 000 in 2 seconds. 



6 



(a)H 



(b) | | 

■ 1 

(c) I I 



(d)-B 




1 



FIG. 6: (color online). The schematic view of the different 
quartz plate configurations. The configurations from (a) to 
(e) correspond to the curves of Fig. 4 of the main text from 
top to bottom. For all cases the maximal length of the quartz 
plates in both arms is 199Ao, giving a maximal total effective 
path difference of 398 Ao- We start to add quartz plates to 
arm 2 when the length of the plate in arm 1 has reached the 
length of (a) 0, (b) 75A , (c) 100 A , (d) 150 A , and (e) 199A . 

Quartz plate configurations 

Figure [6] shows the five different quartz plate configura- 
tions used for Fig. 4 of the main text. To obtain the data 
of Fig. 4, the FWHM of the UV pulse is about 0.18 nm, 
the pump power is about 2.4 mW, and the total coinci- 
dent count rate on the basis vectors HH, HV, VH, VV is 
about 18, 000 in 10 seconds. 
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